home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / clansp.z / clansp
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAANNNNSSSSPPPP((((3333FFFF))))                                                          CCCCLLLLAAAANNNNSSSSPPPP((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLANSP - return the value of the one norm, or the Frobenius norm, or the
  10.      infinity norm, or the element of largest absolute value of a complex
  11.      symmetric matrix A, supplied in packed form
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      REAL FUNCTION CLANSP( NORM, UPLO, N, AP, WORK )
  15.  
  16.          CHARACTER NORM, UPLO
  17.  
  18.          INTEGER   N
  19.  
  20.          REAL      WORK( * )
  21.  
  22.          COMPLEX   AP( * )
  23.  
  24. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  25.      CLANSP  returns the value of the one norm,  or the Frobenius norm, or the
  26.      infinity norm,  or the  element of  largest absolute value  of a complex
  27.      symmetric matrix A,  supplied in packed form.
  28.  
  29.  
  30. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  31.      CLANSP returns the value
  32.  
  33.         CLANSP = ( max(abs(A(i,j))), NORM = 'M' or 'm'
  34.                  (
  35.                  ( norm1(A),         NORM = '1', 'O' or 'o'
  36.                  (
  37.                  ( normI(A),         NORM = 'I' or 'i'
  38.                  (
  39.                  ( normF(A),         NORM = 'F', 'f', 'E' or 'e'
  40.  
  41.      where  norm1  denotes the  one norm of a matrix (maximum column sum),
  42.      normI  denotes the  infinity norm  of a matrix  (maximum row sum) and
  43.      normF  denotes the  Frobenius norm of a matrix (square root of sum of
  44.      squares).  Note that  max(abs(A(i,j)))  is not a  matrix norm.
  45.  
  46.  
  47. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  48.      NORM    (input) CHARACTER*1
  49.              Specifies the value to be returned in CLANSP as described above.
  50.  
  51.      UPLO    (input) CHARACTER*1
  52.              Specifies whether the upper or lower triangular part of the
  53.              symmetric matrix A is supplied.  = 'U':  Upper triangular part of
  54.              A is supplied
  55.              = 'L':  Lower triangular part of A is supplied
  56.  
  57.      N       (input) INTEGER
  58.              The order of the matrix A.  N >= 0.  When N = 0, CLANSP is set to
  59.              zero.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAANNNNSSSSPPPP((((3333FFFF))))                                                          CCCCLLLLAAAANNNNSSSSPPPP((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      AP      (input) COMPLEX array, dimension (N*(N+1)/2)
  75.              The upper or lower triangle of the symmetric matrix A, packed
  76.              columnwise in a linear array.  The j-th column of A is stored in
  77.              the array AP as follows:  if UPLO = 'U', AP(i + (j-1)*j/2) =
  78.              A(i,j) for 1<=i<=j; if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) =
  79.              A(i,j) for j<=i<=n.
  80.  
  81.      WORK    (workspace) REAL array, dimension (LWORK),
  82.              where LWORK >= N when NORM = 'I' or '1' or 'O'; otherwise, WORK
  83.              is not referenced.
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.